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The dynamics associated with the jettisoning of a strap-on booster from the parent launch vehicle is investigated. 
The present analysis accounts for experimental aerodynamics of separation utilizing a two-step method that 
employs a wind tunnel in a problem-solving mode. The first step consists of collecting wind-tunnel test results at 
a large number of points in the multidimensional space of independent variables. A suitable approach for surface 
fitting to the data leads to empirical models of the aerodynamic behavior, which are used in numerical simulations 
to design the appropriate separation system. The suitability of this design is examined, as a second step, through 
extensive simulations using the wind tunnel as a hardware subroutine to generate aerodynamic forces and moments. 
Two strap-on booster configurations are used for illustration. 
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Nomenclature 

= / th aerodynamic coefficient for the yth body, 
with ; from 1 to 6 

= reference values for distance and surface area, 
respectively, for nondimensionalizaticn 
= resultant force and moment components, 
respectively, on the body j along the axes 
= principal mass moments of inertia 
of the y th body 

= indices representing the coordinate axes 
and the bodies, respectively 
= index 

= max(4,4) 

= distance of the base of the core from its own 
mass center; Fig. 1 

= distances of nose of the booster and its base 
from its own mass center 
= freestream Mach number 
= mass of the jth body 
= number of bodies including the core 
= dynamic pressure 

= cylindrical coordinates of the instantaneous 
mass center of the body j with respect to 
C frame 
= tailoff thrust 

= time since initiation of separation 
= velocity of the y th body with respect to the C 
frame expressed in B j frame 
= local core inertial frame, C frame 
= local strap-on inertial frame, S J frame 
= body fixed frame associated with the j th body, 
B J frame 

= instantaneousposition of the mass 
center of the j th body referred to S J frame 
= rig setting parameters; Fig. 3 
= two angles representing the orientation of the 
longitudinal axis of the yth strap-on with 
respect to the C frame; Fig. 1 
= pitch and yaw angle of attack, respectively; 
Fig. 1 
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= rig setting constants; Fig. 3 
= I-evi-Civita density function 
= Eulerian angles representing the attitude 
of the body j 

= rig setting angles; Fig. 3 
= angle of incidence and wind plane angle, 
respectively; Fig. 1 
= spring precompression 
= model scale for scaled down model, > 1 
= angular velocity 


= aerodynamics 
= maximum value of () 

= reference value of () for 
nondimensionalization 
= value of () at the initiation of separation 
= separation system 
= thrust and/or residual thrust 
= change in () 


Introduction 

T HE general dynamics of separating bodies has received the 
attention of several investigators. 1-6 The problem of strap-on 
booster separation, however, has received little attention. Although a 
number of countries have been using strap-on technology for years, 
no directly related investigations appear to have been reported in 
the open literature until recently. Even though Prahlad 7 and Sun- 
daramurthy et al. 8 have touched on some of its aspects, the first 
systematic attempt to analyze this problem was made by the present 
authors. 9 Analyzing the geometry of collisionin a multibody system, 
the authors have developed three different approaches for detection 
of collision 9 and have constructed a no-collision domain in design 
par ameter space. A novel approachis adopted to decouple separation 
dynamics from the core control dynamics, leading to considerable 
simplification. The present paper represents a systematic attempt to 
incorporate the important aspect of separation aerodynamics in the 
design analysis that has previously been ignored. 

The aerodynamics of separation is rather complex due to rela¬ 
tively large angles of incidence, intricate geometry, and the involved 
interference phenomenon. Recently, there have been some attempts 
to simulate this complex flowfield. Moraes et al. 10 developed the 
Navier-Stokes code that provided a good agreement with the ex¬ 
periments but only on the nose part of the on-going stage. A more 
successful attempt seems to have been made by Singh et al., 11 who 
employed a second-orderexplicitpredictor-correctormethod using 
MacCormack's shock-capturing scheme together with overlapping 
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grids to solve Euler’s equations for zero angle of incidence. Evi¬ 
dently, these results cannot be applied in the present situation be¬ 
cause the angles of incidence can be significantly high, further ac¬ 
companied by changing orientations of the separating bodies. This 
as well as similar subsequent attempts at purely theoretical treat¬ 
ment of separation aerodynamics, 12 therefore, appear to be totally 
inadequate. 

A considerable amount of aerodynamic data 13,14 including sta¬ 
bility derivatives 15 has been generated for Space Shuttle-type con¬ 
figurations involving winged bodies. The data have been utilized 
for separation dynamics investigations. 16,17 Lanfranco 18 has inves¬ 
tigated dynamics as well as aerodynamics of separation using a 
captive trajectory system. These earlier investigations are related to 
the separation of two similar-sizedbodiesof fully reusable, parti ally 
reusable (like the current Space Shuttle), and ferry types of config¬ 
urations. These have been effectively summarized by Decker and 
Wilhite. 19 In contrast, not much is available in the open literature on 
aerodynamics of strap-ons separating from the conventional launch 
vehicles. The aerodynamic characteristicshere are far different from 
those of the reusable vehicles. This problem is also different from 
separating an external store from the parent aircraft. Although the 
experimental results provided by Moraes et al. 10 come closer to the 
geometry at hand, the emphasis lies mainly on the aerodynamic 
characteristics of the core in the presence of boosters. This, how¬ 
ever, is of little interest to the present problem. From a separation 
viewpoint, more pertinent test results comprising axial force, side 
force, and yaw moment coefficients for two equal-sized boosters 
have been presented by Sundaramurthy et al. 8 

Formulation 

The details of formulation of the equations of motion are given 
in Ref. 9. Therefore, it is sufficient to state the governing equations 
of motion as follows: 





Fig. 1 Description of various coordinate frames and aerodynamic state 
variables. 


t)/ — e M v{^ = Fj jm’ i = 1, 2, 3 (1) 

ij 6yl - e M ll ccfj co J £ = M] i = 1,2,3 (2) 

where the external force F! and moment M 7 are given by 

A - '•;!' f C *' = i.2,3 

Mj = M 7 + M 7 + M 7 . j = 1, 2, 3,..., n 


Figure 1 describes the various aerodynamic variables required to 
express the aerodynamicforcesandmoments. The externalflowfield 
is specified in terms of two angles, ft, and <P W . Of these, ft, called 
the angle of incidence, represents the included angle between the 
X C i axis and the relative wind velocity vector, whereas the wind 
plane angle <D U , denotes the angle made by the velocity component 
in the X C2 X C3 plane with the X 2 axis. The pitch angle of attack a„, 
and the yaw angle of attack ft, are related to these two angles as 


and models for forces and moments due to thrust, the residual thrust, 
and spring force were developed in Ref. 9. 


Aerodynamics of Separation 

The aerodynamic forces and moments on the body j along the 
t'th axis of its own body frame are expressed as follows: 


f 7 = q w s k{ c; 1 
M 7 . = Q w S Ki D K f C/ +3 J 


j = 1,2,3. n 

i = 1, 2, 3 


(4) 


where D ref and ftef. used as convenient references for distance and 
surface area, respectively, correspond to the strap-on diameter and 
its cross-sectionalarea. The aerodynamic coefficients C 7 denote the 
force coefficients for the values 1,2, and 3 for i , whereas the values 
4, 5, and 6 represent the moment coefficients corresponding to the 
Xi, X 2 , and X 3 axes, respectively (see Appendix). The entire prob¬ 
lem of separation aerodynamics thus reduces to providing models 
or algorithms for computation of these coefficients. 

Several simplifying assumptions are made to render the problem 
tractable. Freestream Mach number, dynamic pressure, and external 
flow are assumed to remain unaltered during the separation process. 
Hence, the change in the angle of incidence results only from the 
change in the orientation of the separating body. A rather short du¬ 
ration of the separation process spanning less than half a second 
justifies these assumptions. Furthermore, the minor asymmetries 
in the geometry of the strap-on due to attachments^trotrusions are 
ignored. 


a u , = arctan(tan ft, cos <P W ) ft, = arctan(tan ft, sin & w ) (5) 


Inversely, 


3 [( 


6 W = arctan (tan 2 a w + tan 2 


)\ 


<p„, 



(6) 


Although it is possible to estimate the bounds on ft, in flight, the 
estimation of the wind plane angle may not be easy as the latter is 
governed by the attitude errors and local wind conditions prevailing 
at that altitude. In addition to these two parameters of the overall 
system, two parametersdescribingthe orientationof the longitudinal 
axis of a booster prove useful in subsequent analysis. Let these be 
denoted as a 7 and ft 7 . These are now expressed in terms of the 
corresponding attitude parameters 



ft = arcsin(cos ft sin ft) (7) 


Here, attention is focused on developing a sufficiently accurate de¬ 
scription of the aerodynamic forces and moments in a suitable form 
facilitating design of the separation system. For this specific objec¬ 
tive at hand, a more general and precise treatment for a complete 
comprehension of the complex aerodynamic issues is not essential. 

For discussion of the methodology, we consider a core with two 
boosters with all of the bodies having the same diameter in a super¬ 
sonic flow. Each of the strap-ons is assumed to be attached to the 
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core through a set of four springs, two near the top end and two near 
bottom. The design problem can now be stated as follows. 


Design Problem 

The design objective is to determine the spring characteristics 
such that for 


M min < < M m 

Umax — Pw — Umax 


Umax — U (/ < 


T < T 
J off — 1 max 


Qw < Qrr 


( 8 ) 


the separation of the strap-ons from the core is ensured without any 
mechanical interference. This paper aims at evolving a strategy for 
design of the springs so as to meet this requirement. 

Evidently, the aerodynamic coefficients for, e.g., the yth strap- 
on, would be influenced by the aerodynamic state variables of other 
strap-ons as well. Hence, in general, symbolically 

Cj = Cj (M K , 9 w , <I> w ,A r J , A , A h J , a J s , 

p^Ar k ,A X k ,Ah k ,a k ,p k ) (9) 

with k representing all of the other strap-ons. Naturally, such a re¬ 
lationship depends on the geometries of the separating bodies. 

To deal with the various issues involved, a two-step procedure 
is proposed. The first one, hereafter called grid approach, consists 
of aerodynamic data generation through wind-tunnel tests by mak¬ 
ing measurements at judiciously chosen discrete aerodynamic states 
and evolving a reasonably accurate empirical model representingthe 
aerodynamicbehaviorof the separating strap-ons. This model, then, 
is utilized to arrive at a suitable design for the springs, through ex¬ 
tensive numerical simulations. This is followed by the second step, 
hereafter referred to as the time march approach, which attempts at 
verification of the design through more accurate simulations using 
the wind tunnel as a hardware subroutine. This technique is known 
as the semicaptive trajectory technique. 20 

Grid Approach 

The grid approach involves building an empirical aerodynamic 
model from the experimentally determined data at well-designed 
grid points of independentaerody namic states. We aim at developing 
the two extreme trajectories within which the actual path of the 
separating strap-ons will lie. 

To limit the wind-tunnel run time to a minimum, several sim¬ 
plifying assumptions are made. All strap-on separation events are 
in supersonic regimes, far away from the transonic Mach numbers. 
Furthermore, the Mach number variation itself is usually small and, 
hence, has little effect on the aerodynamics. 

It may be seen that, for the system under consideration^ strap- 
on may collide with the core in two possible ways, with either the 
nose or the base of the strap-on hitting the core. Therefore, one aims 
at synthesizing the extremes of these two cases. In keeping with 
the earlier approach of a simplified treatment in the first phase, we 
assume the wind velocity vector to be lying in the plane X CI X C2 
as the basis for a conservative design. Therefore, it is apparent that 
aerodynamic out-of-plane forces and moments would be negligible. 
Also, the misalignments of the strap-on tailoff thrust vector and 
that of the separation mechanism would be confined to the yaw 
plane for the synthesis of the most adverse situation. As a result, the 
dynamics under consideration becomes planar. Further, in general, 
the aerodynamics is expected to change monotonically with the 
angle of incidence, as is later confirmed through experiments, and 
therefore, it would suffice to consider only the two extreme angles 
of incidence. Thus, if we can ensure collision-free separationfor the 
case <D W = 0 and 9 W = 0 max that is likely to lead to avoidance of 
collision in virtually all cases involved, then 
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Fig. 2 Influence of the wind plane angle over the aerodynamic coeffi¬ 
cients of the separating strap-ons. 


wind plane angle, whereas that for the leeward strap-on shows some 
improvement when & w is different from 0 deg although not exceed¬ 
ing the error bar. The yaw moment coefficient C 6 , however, displays 
an opposite trend, again within the measurement uncertainty. Nat¬ 
urally, worst dynamics of the windward booster will result from 
& w = —37 deg, whereas the same for the leeward will require 
<P W = 0 deg. Nevertheless, we proceed with the assumptionthat the 
yaw plane aerodynamics will not be affected by the pitch angle of 
incidence for the preliminary design purposes. Naturally, only three 
aerodynamic coefficients, Ci, C 2 , and C 6 , and three aerodynamic 
state variables,/! r, p, : , and A h, are now of significance. 

To minimize the wind-tunnel test program, the physical range of 
interest for the various parameters is chosen. A consideration of the 
geometry of the separating bodies leads to the following range for 
A r and A h: 

A r < ■€. and — A h < (11) 

The estimation of the range of P s is much more involved due to 
a large number of factors, e.g., dynamic pressure, geometry, etc., 
that influence its growth in the time interval of interest dictated by 
the mentioned geometrical inequalities. A trial-and-error approach 
is required to ascertain limits on f } s . 

These ranges and the geometric configurations are adequate for 
the design of a rig setting mechanism to realize these parameters 
in a wind tunnel. A simple and cost-effectivedesign was developed 
by Sundaramurthy et al. 21 A schematic is reproduced here in Fig. 3 
defining the rig setting parameters. The relationship between the 
rig setting parameters x r , y r , z r , 9 r , and Yr and the corresponding 
aerodynamic state variables for the configuration under considera¬ 
tion can be easily obtained geometrically: 

9 r — —a s Yr — Ps y r = — A h/ a — 4 sin Yr 


z r = LHo cos^o — r cos x)l o’ — <5i sin 9 r — 4(1 — cos 9,) 


0 < & w < 7i and -9 nmx < 9 U , < 0 max (10) 


— 4 cos Yr sin 9 r J/ cos 6j. 


( 12 ) 


For cases <P W 0 deg, the out-of-plane aerodynamic force and 
moment components induce a general three-dimensional motion; 
nonetheless, it is not expected to make the collision between the 
separating bodies any more likely. From a few experimental results 
shown in Fig. 2, it is evident that the side force coefficient C 2 for 
the windward boosterremains virtually unaffected by the change in 


x r = —A h/ a — 4(1 — cos 9 r ) + (4 — z r ) sin 6, 

+ 4(1 — cos Yr cos 9 r ) 

A typical set of parametric combinations in the A r-[i s plane chosen 
for the wind-tunnel test runs on the basis of several practicalphy sical 
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constraintsis shown in Fig. 4. These appear in the form of a grid. A 
judicious choice of such a grid pattern is crucial to the development 
of an effective empirical model. The various points in the grid are 
so chosen as to capture the essential characteristics, such as far-field 
and near-field behavior, adequately. 

This entire grid pattern corresponds to one particular value of A h, 
e.g., 0. To include the effects of the variationof the longitudinalgap, 
the data for the same pattern may have to be obtained at different 
values of Ah. This not only increases the number of blowdowns 
considerably but also complicates the surface fitting exercise. As a 
first step, it is now proposed to measure the influence of A h at only 
oneA r-ft combination, e.g., initial, and use the same effect all over 
for preliminary design. Although not strictly valid, we accept these 
simplifying assumptions merely as a first step. 

As yet another simplification, although not truly valid (Fig. 5), 
we assume that the aerodynamic behavior of every strap-on remains 
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Fig. 5 Influence on the aerodynamic characteristics of a strap-on due 
to the variation of A r for the other. 


unaffectedby the changes in the position or orientationof the others. 
This assumption, however, would be relaxed later in the second step 
when the validity of the design would be tested. The empirical model 
under consideration then takes the following simplified form: 

Cj = C/(A r, ft, A h) i = 1, 2, and 6 (13) 


at <P W = 0 and ft, = ±ft lax . 

In view of this assumption, it is now possible to maximize the data 
in each blowdown by instrumenting the two strap-ons and placing 
them so as to cover two different grid points. Considerations of 
symmetry then enable transfer of data collected for one strap-on to 
that of the other as 

Cf(A r,A x, ft)k, = Cf(A r,A X , ft)l-«„ 

i = 1-6 and & w = 0 (14) 

Surface fitting to the ensemble of the aerodynamic coefficients over 
the numerous grid points as already discussed, we obtain the fol¬ 
lowing empirical relations for <P W = 0 deg required for numerical 
simulations. 

Strap-on 1 (body 2): 

C\ = -1.1 -0.02Ar -0.01ft 
C\ = 0.23 + 0.25A r - 0.04A r 1 - 0.002A r 3 
+ 0.125ft + 0.005ft 2 (15) 

C 2 = 1.05 + 2.23A r - 1.64 r 2 + 0.3A r 3 + 0.1ft 
+ 0.005ft 2 -0.08A h 
Strap-on 2 (body 3): 

C 2 = —1.12 + 0.01A r — 0.005ft 
C 2 = -0.26 - 0.02A r - 0.004A r 2 - 0.0007A r 3 

(16) 

+ 0.125ft -0.001ft 2 

C 2 = 0.08 - 1.1A r + 0.5A r 1 - 0.07A r 3 + 0.2ft - 0.05A h 

where ft. is in degrees and all of the distances are normalized with 
respect to D ref . 

The numerical simulations enable us to arrive at a physically 
realizable set of springs, which ensure a collision-free separation 
with a wide margin of safety. Needless to say, the design must ensure 
that the nose, as well as the base of either of the strap-ons,does not hit 
the core even if the uncertainties in aerodynamic measurements and 
the resulting aerodynamic model combine with the environmental 
parameters in the worst possible way. 
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For a system with parameters as shown in Table 1, the applica¬ 
tion of the grid approach results in a realizable separation system 
described in the same table. The scenario of the falling strap-ons is 
as shown in Fig. 6. To show the effect of aerodynamics, the case 
with no aerodynamics is included. 

Time March Approach 

The set of physically realizable spring parameters obtained 
through the grid approach must now be validated, especially in view 
of the various simplifying assumptions made in the first step. The 
initial data together with the system parameters thus obtained serve 
here as the inputs required to initiate an accurate simulation for this 
purpose. The simplifying assumptions are no longer needed. The 
wind-tunnel run for the model configuration corresponding to this 
known initial state generates the aerodynamic data to be utilized 
off-line for the computer run at this initial time step, which in turn 
enables computation of the rig setting parameters corresponding to 
the new configuration of the separating bodies for the subsequent 
blowdown. Thus, the stage is set for repeating this process for the 
next step, which begins with the wind-tunnel run after manually 
changing the model setting to the one just predicted. This process of 


Table 1 Parameter values for the illustrative example 


Configuration: Core with two boosters, all of equal diameters 
Mass of core 17,000 kg 

Thrust on core (constant, corrected 365 kN 

for atmospheric pressure) 


Mass 

Mass moment of inertia 
Leftover thrust 
Nozzle cant 
Thrust location 


Each strap-on booster 

2,900 kg 

700, 33,000,33,000 kgm 2 
3 kN 
10 deg 

4 m below the mass center 


-£{) = 9.35 m £[ = 6.25 m £ 2 = 4.25 m 

7*0 = 1.15 m h 0 = -4.75 m 


Aerodynamics 

Q = 55 kPa S re f = ^/4m 2 D re f = 1 m 

M x = 2.4 Cl = -0.75 

Range of angle of incidence, deg — 3 < a w <3 and —4 < Pw < 4 


Separation system (derived using grid approach and 
subsequently confirmed through time march approach) 
Number of identical springs for each strap-on = four; two near the top 
and two near the base 

Load deflection relation for each spring ( P , kN and 5, m) 

P(S) = 3.5 + 7505-80^ 

0 < 5< p and p = 0.11 m 



a) With aerodynamics (Q = 55 kPa) 


marching in time from an earlier step to the next is continued until 
we are able to ascertain separation with or without collision (Fig. 7). 

The interval for each time march step must be carefully chosen. 
It decides the instants at which the aerodynamic data are sampled 
through wind-tunnel test runs. For the gradually changing system 
configuration, even larger step size may suffice. In case of rapidly 
varying configurations, however, the aerodynamic data must be sam¬ 
pled more frequently and, hence, relatively smaller step sizes must 
be used. Thus, the step size may depend on the initial system pa¬ 
rameters as well as the time elapsed during the separation process. 
Moreover, smaller step size, although capable of capturing aerody¬ 
namic behavior more faithfully, would be more susceptible to errors 
in model setting and the resulting measurement noise. A suitable 
compromise needs to be struck between these conflicting require¬ 
ments. For the systems considered here, initially, a step size of 40 
ms was found suitable. For reasons explained earlier, this had to be 



Fig. 7 Flow chart showing sequence of steps in time march approach. 




b) Without aerodynamics (Q = 0) 


Fig. 6 Trajectories for the separating strap-ons. 
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gradually reduced to 20 ms at the end of the simulation. The stability, 
convergence, and accuracy of the numerical integration of the gov¬ 
erning equations require the integration step to be much smaller (^ 1 
ms) compared to the step size at which blowdowns are conducted. 
Therefore, some suitable interpolation/extrapolation schemes would 
also be needed to compute required aerodynamic coefficients at de¬ 
sired close integration steps from blowdowns conducted at coarse 
time march steps. A judicious choice of reference variable for this 
purpose enables a more effective utilization of the sampled data. 
The comparisonbetween the extrapolationand the subsequentblow- 
downprovidesan in-built mechanism to regulate the time march step 
size. 

The numerical simulation duly accountsfor the core thrust and the 
tailoff thrust of the strap-ons in addition to the aerodynamic effects 
discussed earlier. The pertinent details of the models for separation 



a) <F„ = 0 deg b) = —37 deg 

Fig. 8 Simulation results for a two-booster configuration confirming 
collision-free separation. 


systems, residual thrusts, and the methods of detection of collisions, 
if any, used in the analysis are available elsewhere. 9 

Results and Discussion 

For the results presented here, the wind-tunnel tests were con¬ 
ducted at the 1.2-m trisonic facility of National Aerospace Labora¬ 
tories. Bangalore, India. For the system under consideration, the sys¬ 
tem parameters are defined in Table 1, wherein the separation system 
obtained from the grid approach is also included. The simulation re¬ 
sults obtained with a view to validating the configuration developed 
in the first step are generated through a time march approach. The 
results are provided for two wind plane angles (Fig. 8). If desired, 
more cases can be examined to achieve full confidence in the design. 

The simulations undertaken here clearly suggest ever present 
large clearances that ensure a collision-free separation with a wide 
margin of safety. The angle <P W = 0 deg corresponds to the case of 
relative motion in a single plane. However, the case <P W = — 37 deg 
(Fig. 8b) leading to a general three-dimensional motion demands 
the clearances in both the views as shown here. It may not always 
be possible to ensure a nose opening rotation for the boosters under 
all possible circumstances due to various geometrical and physi¬ 
cal constraints on the implementable separation system. Although 
a nose closing rotation of the windward booster is evident here, a 
minimum physical clearance, under the worst possible combination 
of environmentalparameters and the aerodynamic modeling uncer¬ 
tainty, beyond 0.25Z) re f was considered adequate. The presence of 
additional attachments on the core, e.g., fins, enlarges the prohibited 
region that must not be violated by the strap-on axis of symmetry. 
This enlargement can be obtained through simple geometric con¬ 
sideration of the sizes and the core attitude excursion, if any. 

On the other hand, should any simulation run lead to the occur¬ 
rence of collision between any two separating bodies, the prelimi¬ 
nary design of the separation system should be suitably altered and 
subsequently confirmed using the described approach. 

In the second example, we considera more complex situation of a 
launch vehicle with six strap-ons, of which only four are jettisoned. 
Unlike in the preceding configuration, here we must also ensure 
that the separating strap-ons do not hit the nonseparating adjacent 
ones. The number of modes of possible collision has now gone up 
to four. Accordingly, while using the grid approach of the first step, 
additional wind plane angles need to be considered to cover the new 
possibilities. Figure 9 shows the results of the simulation using the 
time march approach for one wind plane angle only. Collision-free 
separation is evident for this case. Simulation runs for other criti¬ 
cal cases, not shown here, also confirm the suitability of the design 
chosen. 

The suitability of the preliminary design obtained in the first step 
has been invariably confirmed. The methodologyproposedhere has 
already been proven in successful flights of the Indian launchers 



Fig. 9 Simulation results for a six-booster configuration confirming collision-free separation. 
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ASLV (Augmented Satellite Launch Vehicle) and PSLV (Polar Sat¬ 
ellite Launch Vehicle). 


Appendix: Frames of Reference 

The coordinate systems used for convenience of representation 
and interpretation are described. 

1) Local core inertial frame X CI X C2 X C3 (C frame) with origin 
frozenatthecenterofmassatt = 0 and taking the axis X C1 alongthe 
longitudinal axis, X C1 along the pitch axis, and X C3 was chosen so 
as to complete the right-handedtriad (Fig. 1). This frame is assumed 
to move with uniform system velocity taken at t — 0. Because, at 
this instant, the core and the strap-ons have the same orientation 
and velocity with reference to this C frame, the initial conditions 
of velocity and orientation for all of the bodies reduce to zero. The 
velocity of the core (and all other bodies), however, may evolve with 
time under the influence of external forces. 

2) Local strap-on inertial frames X J S1 X J S2 X J S3 (S J frame) for the 
j th strap-on, with origin at its own center of mass at t = 0, was 
takenparalleltoX cl , X J S2 along the radial direction,and X J S3 was se¬ 
lected so as to complete the right-handedtriad as before (Fig. 1). 

3) Body frame X[X 2 X 3 (B j frame) for the jth body was such 
that the frames B j and S J are identical at t = 0. 

The body frames are introduced to facilitate the Eulerian descrip¬ 
tion of angular motion. These noninertial frames are obtained from 
the inertial S J frames (C 1 frame) by three successive rotations in the 
following sequence: 1) a rotation of 0 3 about the X J S3 axis (Xc 3 axis) 
leading to intermediate frame (X' X' 2 X' 3 ), 2) a rotation of 6^ about 
the X' 2 axis leading to anotherintermediateframe (X ^ X 1 X 3 ), and, fi¬ 
nally, 3) a rotation of &[ about the X t axis leading to the body frame 

The origin of the S J frame, the mass center of the body j at 
t = 0, is referred to the C frame in terms of the radial distance r () ', 
the angular distance Xg, an d the longitudinal location h J 0 . The S J 
frame is obtained from the C frame by a single rotation of Xg about 
the X cl axis leading to a time-independent transformation matrix 
relating the two. The parameters, representing the mass center of 
the body j at some general time t, are referred to as r' , , and h '. 

The instantaneous change in the position of the mass center of the 
jth body, represented by A r\A x J , and A h', is then given by 

A > J = {7(''o cos ^o + x i) 2 + ( r o sin ^o +-* 3 )’} ~ r o 


A x 1 = arctan 


rj sin Xp + *3 | 
Jg COSXg +X J 2 \ 


■ Xo 


A h,i = x[ — h J Q 
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